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ABSTRACT 

The low energy spectral slopes of the prompt emission of most gamma-ray bursts (GRBs) are difficult to rec- 
oncile with radiatively efficient optically thin emission models irrespective of the radiation mechanism. An 
alternative is to ascribe the radiation around the spectral peak to a thermalization process occurring well inside 
the Thompson photosphere. This quasi-thermal spectrum can evolve into the observed non-thermal shape by 
additional energy release at moderate to small Thomson optical depths, which can readily give rise to the hard 
spectral tail. The position of the spectral peak is determined by the temperature and Lorentz factor of the flow 
in the termalization zone, where the total number of photons carried by the jet is established. To reach ther- 
malization, dissipation alone is not sufficient and photon generation requires an efficient emission/absorption 
process in addition to scattering. We perform a systematic study of all relevant photon production mechanisms 
searching for possible conditions in which thermalization can take place. We find that a significant fraction 
of the available energy should be dissipated at intermediate radii, ~ 10^°- a few x 10^^ cm and the flow there 
should be relatively slow: the bulk Lorentz factor could not exceed a few tens. The least restrictive constraint 
for successful thermalization, F < 20, is obtained if synchrotron emission acts as the photon source. This 
requires, however, a non-thermal acceleration deep below the Thomson photosphere transferring a significant 
fraction of the flow energy to relativistic electrons with Lorentz factors between 10 and 100. Other processes 
require bulk flow Lorentz factors of order of a few. We examine the implications of these results to different 
GRB photospheric emission models. 

Subject headings: gamma-ray burst: general — gamma rays: general — radiation mechanisms: thermal — 
radiation mechanisms: non-thermal — radiative transfer — scattering 



1. INTRODUCTION 

The typical broadband spectra of gamma-ray burst (GRB) 
prompt emission can be broadly characterized by the follow- 
ing properties: a peak at a few 100 keV, a non-thermal power- 
law extending from the peak to higher (in some cases GeV) 
energies, and a rising but non-thermal slope below the peak. 

The enormous luminosity coupled with the non-thermal ap- 
pearance of the spectrum naturally suggests that the emis- 
sion is produced by non-thermal high-energy particles ra- 
diating in an optically thin environment. However, opti- 
cally thin emission models face serious difficulties in explain- 
ing the low-energy turnover i n the obse rved spe ctra, regard- 
less of the emiss ion process jCohen etal. 1997: .Crider et al.l 
II997HPreece"era l. 1998; Ghir landa et ani2003l) . This has led 
to models in which energy dissipation takes place deep in- 
side the Thomson photosphere, where the thermalization pro- 
cesses are efficient enough to produce a close to blackbody 
spectrum jEichler & Le vinson 2000; Meszaros & Rees 20()ft 
Rees & Mes"zarosll2005l; [Thompson et alJi2007l; IBeloborodovl 



2O10ir In such models, the radiation field, which has a quasi- 



thermal appearance well below the photosphere, is distorted 
into a non-thermal shape by additional heating and Compton 
scattering near the optically thick-thin transition. 

An important question in photospheric GRB emission mod- 
els concerns the origin and amount of photons present in the 
flow at the stage when it approaches transparency. In the sim- 
ple fireball model these photons are created near the central 
engine and advected with the jet towards the photosphere. The 
typical energy available per photon at the launching site of the 



jet is ^ 5MeV. I n a pure pair-photon fireball jGoodmanl 19861 ; 
|Paczvnskill986l) this would also be the average observed pho- 
ton energy. On the other hand, a significant baryon load can 
degrade the mean photon ener gy by converting mos t of the 
internal energy to bulk motion TShemi & Piranlll990l) . How- 
ever, if no more photons are created in the jet along its way, 
any process making the flow radiatively efficient would bring 
the average photon energy back to a few MeV. The disagree- 
ment with the observed spectral peak positions thus requires 
an additional dissipation/photon production far from the cen- 
ter. 

The frequently used assumption that any dissipation tak- 
ing place well below the Thomson photosphere leads to 
a thermal (Planck-like) radiation spectrum is generally in- 
valid. The crucial ingredient is the existence of an emis- 
sion/absorption process allowing the radiation field to attain 
a therm odynamic (rather th an just kinetic) equilibrium with 
matter (IBeloborodovl 120 1 2h . Without the creation of addi- 
tional photons (as in Compton scattering) sub-photospheric 
dissipation can only redistribute energy between bulk mo- 
tion/Poynting flux and the existing photons, while the total 
flow energy per photon remains unchanged. Unless the ra- 
diative efficiency is low, the spectra peak at too high ener- 
gies compared with observations. The peak energy can only 
be decreased by introducing new photons. The lowest attain- 
able energy is set by the blackbody limit reached upon com- 
plete thermalization (E ichler 1994; Eichler & Levinson 200(| 
iThompson etani2007l) . 

The purpose of this work is to study the efficiency of photon 
production and thermalization deep inside the photosphere of 
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GRB outflows as well as the resulting spectral peak energies, 
accounting for all plausible emission/absorption processes. 
For GRB jets the possible candidates are cyclo-synchrotron 
emission, bremsstrahlung and double Compton scattering. 
Subject to the existence of a photon source, the relaxation of 
the radiation field to a Planck spectrum can be achieved in two 
ways: the most straightforward option is for the jet to be opti- 
cally thick to absorption throughout the spectrum. The second 
and less restrictive option can be realized if the photons pro- 
vided by the emission process are redistributed into a thermal 
spectrum by saturated Comptonization. In the latter case two 
conditions have to be satisfied: the number of produced pho- 
tons must be sufficient to fill the Planck spectrum, and the 
Compton y-parameter has to be large for Comptonization to 
proceed in a saturated regime. 

We examine the problem both analytically and numeri- 
cally. For the latt e r we use the kin etic code develope d by 
IVurm & PoutanenI (120091) (see also IVurm et alJ 1201 Ih that 
self-consistently solves the coupled kinetic equations for elec- 
trons and photons in relativistic flows. For the present prob- 
lem the code was modified to include induced Compton scat- 
tering, which was missing in the original version. 

The paper is organized as follows: We begin by briefly out- 
lining the spectral hardness problem of optically thin emission 
models. We continue with a general discussion of photon pro- 
duction and thermalization in the jet as well as their relation 
to the observed spectral peak energies in the context of opti- 
cally thick emission models. This is followed by a detailed 
study of various processes acting as photon sources in GRB 
jets. Finally, we discuss the results and their imphcations to 
GRB models. 

2. THE LOW-ENERGY TURNOVER PROBLEM 

The prompt emission from GRBs is observed in a wide fre- 
quency band, from hard X-rays to hard gamma-rays. The 
energy spectrum, E^dN^/dE, peaks at Epi_ ~ 100 — 
1000 k eV. The ob served sp ectral peak-luminosity rela- 
tion is dYonetoku et al. 2004; iGhirlanda et al] l2009l |20T1 
iGruber et al.ll20Tll) 

Ep^ = 300Llil,,keY, (1) 

where Liad is the observed isotropic equivalent luminosity 
and the notation A = lO'^Ax in cgs units has been used. 
Above the peak, the hard tail with the photon index a bit 
steeper than a = —2 extends till hundreds of MeV (and even 
a few GeV in some cases), which implies a non-thermal emis- 
sion mechanism, e.g. synchrotron or inverse Compton. 

The distribution of ph oton indices below t he peak is cen- 
tered at a w -1 (e.g. iKaneko et aTl 120061: [Goldstein et al.l 
12012*) and extends to values as hard as a = 1 or even beyond 
^hirlanda et al. 2003). The hard low-energy spectra pose a 
great challenge to models in which the whole spectrum is pro- 
duced by relativistic electrons radiating in the o ptically thin 
regime. The synch rotron line of death, a > —2 /3 dKatzl 19941: 
ICohen et al.lll997h, is violated by significant fraction of bursts 
jPreece et al.ll 19981) . which thus defy explanation by optically 
thin synchrotron emission. Furthermore, high radiative effi- 
ciency in the prompt phase implies that electrons emit in the 
fast cooling regime. Independently of the emission mecha- 
nism, this leads to strong excess in the low-energy part of the 
spectrum, in apparent contradiction with observations. This is 
the low-energy turnover problem. 

The problem can be understood in general terms by consid- 



ering any radiation mechanism in which the typical energy, 
E, of the emitted photon scales with the electron's Lorentz 
factor, 7, as i? oc 7''. An electron cooling by A7 emits 
the energy /e^E — mcC^A-f, where /e is the fluence per 
photon energy interval. The fluence due to a single elec- 
tron cooling all the way down from its initial energy is then 
/e = rriec'^d'y/dE oc E^^^^^. Since the cooling time is 
short, all accelerated electrons will cool. The emissivity at 
a given energy E is proportional to the rate at which elec- 
trons are injected above the Lorentz factor 7 corresponding 
to E, which equals to the flux of electrons through 7 towards 
lower energies. If the injection rate is N-mj 7"'', this flux 
is oc 7^P+i oc if p > 1, and it is constant for 

harder injection spectra {p < 1). Correspondingly, the emis- 
sivity is Je oc /eE'^-p+'^^/p oc £;(2-p-rt/p or Je oc E^/p-^ 
in the cases p > 1 and p < 1, respectively. Note that except 
for the requirement of fast cooling, the above argument is in- 
dependent of the interaction cross-section or electron cooling 
rate. Assuming hard electron injection (p < 1), the hard- 
est spectra that can be expected from synchrotron or Comp- 
ton scattering in the Thomson regime (p ~ 2) have a photon 
index a = 1/p — 2 = —1.5, whereas for Compton in the 
Klein-Nishina regime and bremsstrahlung (p = 1) one finds 
a = — 1. Therefore spectra harder than a > —1 are very 
difficult to produce by optically thin emissiorfl 

3. THE PHOTON PRODUCTION PROBLEM 

The low-energy turnover problem is avoided if the low- 
energy spectrum is shaped in the optically thick regions of 
the flow. Such photospheric emission models usually invoke 
some form of dissipation close to the Thomson photosphere 
to account for the non-ther mal app earance of the ov erall spec- 
trum (e.g.'Rees & Meszaros'2005;'Pe'er et al.'2006VGiannios| 
2006 ; Giannios & Spruit 2007; Giannios 2008; Beloborodo^ 
2010l:IVurm et alJI2011l) . However, as shown below, the total 
number of photons carried by the jet is established at much 
smaller radii. Along with the luminosity this number deter- 
mines the position of the spectral peak, compatibility with the 
observed i?pk values thus requires sufficiently efficient pho- 
ton production in the flow. 

For a given luminosity, radius and Lorentz factor the lowest 
-Epk value is attained if the spectrum is a blackbody. Consider 
the case where a fraction eeb of the total available energy L 
is dissipated and processed into a Planck spectrum at radius r, 
which we call the thermalization radius. The isotropic equiv- 
alent radiation luminosity at that location is given by 

eBBL = A^cr^V^^aT^, (2) 

where F and T are are the jet Lorentz factor and (comoving) 
temperature at r, respectively. 

Between the thermalization radius and radius where the ra- 
diation decouples from the flow (i.e. the Thomson photo- 
sphere), the radiation luminosity can be altered by adiabatic 
cooling as well as further dissipation. Let's define eiad as 
the final radiative efficiency, i.e. the fraction of the total en- 
ergy L released to the observer (Lrad = ^rad^)- Without 

' Somewhat hai'der spectra can be produced if different mechanisms are re- 
sponsibl e for electron coolin g and its emission in the spectral range of interest 
(see e.g. IDaigne et alj|201 Ih . A specific mechanism that can keep accelerat- 
ing the same sub-population of electrons can also lead to harder spectra: in 
this case the limit a = —2/3 could in principle be achieved for synchrotron 
emission. 
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further dissipation beyond the thermalization zone, we have 
Eiad < £bb, where the equality holds if the flow remains 
radiation-dominated until the Thomson photosphere. The in- 
equality holds in a coasting flow, in which case the average 
photon energy is degraded by adiabatic cooling. With further 
dissipation one may also have Sj-ad > ^bb- 

If no more photons are produced beyond the thermalization 
zone, the observed spectral peak is approximately given by 



E'pk = GkBTT 



Erad 



(3) 



where GfeeTT is the approximate vFi, peak energy 
of a Lorenz boosted blackbody spectrun j^ (see e.g. 
iGrimsrud & WassermanI [T998t iLi & Sarill2008h and the fac- 
tor erad/^BB accounts for adiabatic cooling and/or additional 
dissipation between the thermalization zone and the Thomson 
photosphere. From equations (|2]i and (O we obtain 



= 200,/^^ ^ (eBBi52)'/' keV. 
V ri2 Ebb 



(4) 



Comparing with the observed position of the spectral peak ([T]i 
and using Liad = Erad^, one finds a relation between the 
thermalization ra dius and the flow Lorenz factor in the ther- 
maliz ation zone jEichler & Levinso n 2000; Thompson et al.l 
l2007l) : 



r 



-rad r-1/2 

372^52' cm. 

-BB 



(5) 



An immediate consequence of Equation ^ is that the ob- 
served positions of the spectral peaks cannot be explained by 
relic photons produced near the central engine, unless the ra- 
diative efficiency is very low. Instead, most of the observed 
photons have to be produced far from the center. 

The derivation of Equation Q assumes that the flow main- 
tains a constant opening angle between the thermalization lo- 
cation and the Thomson photosphere. If the flow undergoes 
substantial collimation between these regions, the constraint 
on r/r is somewhat relaxe d (Equation ( IC3b in Appendix ICl 
see also lBeloborodovll20 1 2h . However, for realistic jet open- 
ing angles the constraint on r is still incompatible with the 
size of the central engine, thus our conclusion that bulk of the 
observed photons have to be produced further in the jet still 
holds. 

Collimation that operates after the thermalization zone 
would somewhat relax condition (|5]l. However, most of the 
jet collimation is likely to take place below the thermalization 
zone. We will therefore use Equation Q as it stands and ig- 
nore the effects of possible further collimation at larger radii, 
which could influence the constraint on r/T at most by a fac- 
tor of a few. 

There are several aspects to be considered regarding the 
way thermalization can be achieved in GRB jets. First, one 
has to stress that thermalization requires that blackbody ra- 
diation is formed. For that the thermalization zone should 
be optically thick to absorption. Dissipation alone below the 
Thomson photosphere is insufficient. The photon production 
problem arises because the rates of emission/absorption pro- 
cesses rapidly decrease with the distance and with the Lorentz 
factor of the flow (because all relevant parameters like the 
proper density and the magnetic field decrease). Therefore, 

^ The uFu peak of a blackbody spectrum in the comoving frame is at 
3.92fcBT, which is boosted on average by l.SF in a relativistic quasi- 
spherical flow, yielding a peak at S.SSfce^r in the external frame. 



it is difficult to find a powerful enough photon source and sat- 
isfy Equation ^ simultaneously. This combination requires 
rather specific conditions within in the jet, as we show below. 

Secondly, thermal-like spectra can be formed via Comp- 
toniz ation of soft photons on thermal electrons (e.g. 
i Liang et al.llI997l; lGhisel lini & Celotti 1999; Thompson et al.l 
l2007HGianniosl l2012; Be loborodovfl2012b . In fliis case, flie 
spectrum also peaks at the energy (O, however the energy 
density may be less than aT'^, therefore Equation (|5]i be- 
comes an inequality so that the production of observed pho- 
tons should take place at an even larger distance. On the other 
hand, if the emission/absorption process is able to provide 
enough photons, complete thermalization may be achieved 
and Equation (|5]l remains as it stands. In general, the condi- 
tions are less restrictive for thermalization by Comptonizing 
soft photons than by the absorption process alone. 

The Comptonization efficiency also places limits on the pa- 
rameters in the thermalization zone. The condition that pho- 
tons are efficiently redistributed towards the thermal peak is 



y 



kv.T 



TOcC 



fTTiVctdyn > 10, 



(6) 



where y is the Compton parameter, N is the proper elec- 
tron number density and tdyn = r/cT the proper propaga- 
tion/dynamical time. One can conveniently normalize the 
proper plasma density by the Lorentz factor corresponding to 
the total energy conversion into the kinetic energy, Fmax: 



(7) 



Then the Thomson o ptical depth of the flow is (e.g 
lAbramowicz et al.|[T99l1) 



TT = CTT-'Vcidyn = = 1-2 • 10 



r noF^r 

max, J 



(8) 



The plasma temperature can be expressed via the peak energy 
(l3|. Using the observed E-py - Liad relation ([T]), the condition 
Q can now be written as 



100- 



.3/2 
^52 



4/d-ior?r 



> 1. 



(9) 



max, 3 



Taking into account condition (|5]), this constrains the Lorentz 
factor in the thermalization zone to a few tens: 



5/8 jl/2 
r ^ 91 ^BB ^52 
^ ^ ^\3/8pl/4 ■ 
rad max, 3 



(10) 



Thus, the kinetic energy should be small in that location. We 
note in passing that this could suggest that the energy is trans- 
ferred from the central engine to the thermalization zone in 
the form of the Poynting flux which accelerates much slower 
than a thermal fireball. Alternatively, one has to find a highly 
efficient mechanism to convert almost all of the available ki- 
netic energy to internal form at r ~ 10^° - 10^^ cm. 

Even stronger restrictions are obtained from the considera- 
tion of the soft photon sources necessary for the Comptoniza- 
tion. The condition (|9]l only ensures that the photons are re- 
distributed towards the observed peak energy. However, one 
has to find a photon source capable of producing the observed 
amount of photons. The condition for the complete thermal- 
ization is written as 



tdyn^>A/'l 



BB: 



(11) 
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where A^^ is the photon injection rate per unit volume, and 



BB 



(12) 



is the total number density of photons in the Planck spectrum. 
Here B = k-oT/m^c? is the dimensionless temperature and 
Ac = h/rricC is the Compton wavelength. We show below 
that the requirement that this condition is satisfied simultane- 
ously with Equation ^ places severe limits on the parameters 
of the thermalization zone. 

In unmagnetized flows the only plausible emis- 
sion/absorption mechanisms are bremsstrahlung and 
double-Compton scattering. Their corresponding absorption 
photospheres are deep below the Thomson photosphere, at 
tt > 10"'. Due to Comptonization, thermalization can occur 
above the absorption photosphere but not far from it. In 
magnetized flows, one has also to take into account cyclotron 
and synchrotron emission, which are copious sources of low 
energy photons. In this case saturated Comptonization of soft 
photons is crucial for thermalization. 

4. THE PHYSICAL MODEL 

To Study the thermalization and photon production capabili- 
ties of the jet by various processes, we consider an energy dis- 
sipation episode taking place between 7'i„in and r,„ax = 2rmin 
in magnetized or unmagnetized outflows carrying baryons, 
electrons and radiation (and at some conditions, electron- 
positron pairs). The condition (|5]l requires that a sizable frac- 
tion of the total available energy - kinetic or magnetic - be 
converted to internal (thermal) energy at radii r > 10^" cm. 
Without specifying the conversion mechanism, we assume 
that this energy is supplied to the electrons either in the form 
of heating the bulk of the lepton population, or by accelerating 
a fraction of them to relativistic energies. The electrons trans- 
fer their energy to radiation via emission of photons by vari- 
ous processes and their Comptonization. The main questions 
we ask for each emission process are whether there is time to 
emit and reprocess enough photons to fill the Planck spectrum 
and whether the required conditions for this are compatible 
with the general constraint (|5]l. We also check the possibility 
of partial thermalization when a Wien rather than the Planck 
spectrum is formed with the peak energy compatible with the 
observed relation 

We study the problem analytically as well as numerically in 
more complicated cases involving magnetized jets. The ana- 
lytic study is based on analysis of the Kompaneets equation, 
which allows us to estimate the amount of photons available 
for thermalization without solving the full equation (see Ap- 
pendix |A]i. For numerical simulati ons we use the kin e tic ra- 
diat ive transfer code d eveloped by IVurm & PoutanerJ (l2009h 
and lVurmet"an(l2011h (Appendix |B]i. 

5. THERMALIZATION BY BREMSSTRAHLUNG 

Bremsstrahlung can serve as a photon source if the density 
is sufficiently large. This can happen only if energy is dissi- 
pated deep in the flow where the temperature is high enough 
to sustain a significant population of electron-positron pairs. 
If the bremsstrahlung optical depth is above unity all the way 
to the thermal peak energy E ^ SfceT, complete thermal- 
ization is obviously achieved. However, a somewhat weaker 
condition for thermalization can be obtained by taking into 
account Comptonization of bremsstrahlung phot ons. 

Making use of Equation ( IA13b in Appendix lA.ll for the 



bremsstrahlung photon injection rate, one can write the ther- 
malization condition (fTTT i as 



Mr ^dyn __ 2^ 2.35fcBr ^ ^ 

ATbb ~ 4.8 £^0 ^ ' 



(13) 



where 



31/2 27r3/2 



-7/2 



A£A^_TT (14) 



is the free-free optical depth at E' = fee 2^- The transition 
energy Eq above which bremsstrahlung photons ar e Com p- 
tonized rather than reabsorbed is given by Equation ( IA12I ). 

Condition (fTsT l can be satisfied only at the stage where the 
flow is still heavily loaded with pairs, which are in (quasi- 
)thermodynamic equilibrium with the radiation. If the radi- 
ation field is a blackbody, the electron/positron chemical po- 
nicC^ and the density of electrons/positrons 



tential is ii± 
is 



(15) 



Substituting this together with Equation (fT4l i into Equation 
( fT3] l yields a constraint on temperature where thermalization 
is still effective, in terms of the flow's Lorentz factor and ra- 
dius. For Planckian radiation we can use the relation ^ to 
eliminate either T 01 r and obtain a constraint for the other 
variable in terms of F, L and £bb- 

We find that Comptonization can extend the thermalization 
region up to Tg^±{E = Sk^T) ^ 0.005. However the ther- 
malization radius is only about 50% larger than the radius of 
the bremsstrahlung photosphere, owing to the strong depen- 
dence of the opacity on the temperature and thus also on r via 
Equation Therefore Comptonization is relevant only in a 
narrow range of radii and thus has a limited role in forming 
the Planckian spectrum. 
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Figure 1. Constraints on r and F from the requirement of efficient 
bremsstrahlung thermalization (solid lines), and from the observed peak en- 
ergies (dashed lines), for flow luminosities L = 10^^ (blue), 10^^ (red), 
IqM gj.g g-i (gfeen) and ebb = 1- The allowed region is below the solid 
and above the dashed lines. The bremsstrahlung photosphere is shown by 
dotted lines. 

Figure [T] depicts the allowed parameter region for r and F, 
for fixed L and ebb- The solid lines correspond to the ther- 
malization condition ( fTSl l and limit the dissipation radius and 
Lorentz factor from above. The constraints from the observed 
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peak energy from Equation ^ are shown by dashed lines, 
limiting r from below. The allowed parameter space forms a 
wedge-shaped region to the left of the crossing point of the 
two constraints. This point sets the maximal Lorenz factor 
in the thermalization zone for a given luminosity. One can 
see that a solution exists (i.e. the maximal F > 1) only for 
luminosities above a fewx 10^^ erg s^^. 

The lower limits on the temperature in the thermalization 
zone are shown in Figure |2] The temperature at which ther- 
malization fails depends weakly on other parameters and is 
around 28 keV. This is a consequence of the exponential de- 
pendence of the number of pairs on 6 (Equation (fTSll). An ap- 
proximate analytic thermalization constraint on r and F can 
be read from Equation (|2]i, keeping T constant at 28 keV: 



rrg < 15(eBBi5 



^l/2 



(16) 



This constraint is compatible with the condition (|5]l only if 

(17) 



P < r, Ebb ^1/2 
1 ^ — ^52 ' 



rad 

i.e. the flow can only be mildly relativistic at the thermaliza- 
tion zone. 
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Figure 2. Bremsstrahlung thermalization constraint on the comoving tem- 
perature (lower lines), for flow luminosities L = 10^^ (dashed), 10^^ (solid), 
10^'' erg (dot-dashed) and sbb = 1- The allowed regions for different 
luminosities are above the corresponding lines. Also shown is the tempera- 
ture at the bremstrahlung photosphere (upper lines) for the same luminosities. 



5.1. The effect of clumping 

In a smooth flow bremsstrahlung rapidly becomes insignifi- 
cant as a source of photons once the pairs drop out. However, 
by virtue of the quadratic dependence on density the average 
emission/absorption rates can be increased if regions of the 
flow are locally compressed to higher densities. Such situ- 
ation can arise, for example, in a Poynting dominated flow 
carrying an alternating magnetic field (striped wind). This 
implies the existence of current sheets where matter density is 
substantially enhanced compared to the average in the flow. 

Let's assume that a fraction A of the plasma is compressed 
in the radial direction by a factor ^ = l/A so that A is the av- 
erage size of the compressed region and I is the average dis- 
tance between such regions. The density in the compressed 
regions is A^dcnsc = A^A^, where N is the average density 
given by Equation Bremsstrahlung opacity is propor- 



tional to the emission measure 



<n^ > r/F 



2 A r 2 2^ 
"dense yp = A £,N — , 



(18) 



and is given by Equations (IA8b and (lAlOl ). Radiation is ther- 
mal as long as afftdyn > 1 at E' = Sk^T, yielding 



Frio < 1-7 X 10~3(A2e)^/^ 



L 



9/10 
52 



-7/10 1 
•BB ' 



,8/5 ' 
max, 3 



(19) 



where we have used Equations (|2]i, ^ and (H) to express fceT, 
N and tt, respectively. We see that high compression factors 
are required for compatibility with condition ©. 

In the striped wind, the maximal expected compression can 
be estimated by balancing the magnetic pressure outside the 
current sheets with the thermal pressure inside, B^/Stt = 
2NdonsckBT. Defining as the fraction of total energy car- 
ried by the Poynting flux by 



sbL 



we obtain 

A^den 



TV 



= 4.3 X 10^ 



1/2 

£S f 10 i max, 3 
ry'(£BBi52)l/4' 



Using Equation ( fT9] l, this gives 

X^/' e'J' Lli' 



Tr^!J < 120 ■ 



^52 

9/14 p4/7 
^BB ^ max, 3 



(20) 



(21) 



(22) 



Compatibility with condition (|5]l limits the Lorentz factor to 



F < 35 



\ 2/5 ^2/5 rl/2 



^52 



3/10 p2/5 
rad max, 3 



(23) 



Note that this condition was obtained neglecting Comptoniza- 
tion because the condition for the saturated Comptonization 
dlOl l is stronger 

6. THERMALIZATION BY DOUBLE COMPTON SCATTERING 

Double Compton scattering, where an additional photon is 
produced upon the scattering of an electron and photon, can 
supply enough photons for thermalization if the dissipation 
operates deep enough below the Thomson photosphere. Mak- 
ing use of th e corresponding photon injection rate ( lAlSI l (Ap- 
pendix|Aj2]), one writes the thermalization condition (fTTl i as 



A/dc idyn 16af 



BB 



'.9Dc(e) In 



kBT 
Eo 



> 1, 



(24) 



where the transition energy Eq is given by Equation ( IA17I ). 
The blackbody relation (|2J| can be used together with Equation 
^ for Tt to eliminate one of the variables 6*, r or F from 
the thermalization condition ( l24l i. Eliminating either r or F, 
we find that temperature at the thermalization limit is a weak 
function of the other variables, scaling approximately as 



9. oc L-V8ri/4ri/4 



-l/6pl/6 



(25) 



and has value around 5 - 10 keV for typical parameters (see 
also lBeloborodovll20 1 2h . Consequently, Equation ( IA17I ) im- 
plies that the logarithmic factor in Equation ( l24l i is approxi- 
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mately constant, lii(fcBr/i?o) 
comes 



4. The constraint itself be- 



rrJ/3<8.5^ 



1/6^1/2 



52 



r. 



/3 



max, 3 

Together with condition (|5]l this requires 



r < 4.1 



7/10^1/2 
^BB 

2/5pl/5 ■ 
rad max, 3 



(26) 



(27) 



Figure [3] depicts the allowed parameter space r and F 
for double Compton. The thermalization constraint ( |26] | is 
shown by solid lines and the dashed lines correspond to 
the constraint from the observed peak energies. Equation 
Q. Comparison with Figure [T]reveals that thermalization is 
achieved via the double Compton scattering a bit easier than 
via bremsstrahlung. Still the conditions are rather stringent; 
the process works efficiently only if the energy is released in 
a mildly relativistic flow. 




Figure 3. Constraints on r and T from the requirement of efficient thermal- 
ization by double Compton scattering (solid lines), and from observed peak 
energies (dashed lines). The blue, red and green lines correspond to luminosi- 
ties L = 10^^, L = 10^^ and L = 10^"^ erg s^^, respectively. In all cases 
£bb = 1- For Each luminosity the allowed region is below the corresponding 
solid and above the corresponding dashed line. 



7. THERMALIZATION BY COMPTONIZING CYCLOTRON 
RADIATION 

In highly magnetized flows, one has to take into account 
cyclotron radiation of thermal electrons. If the fraction £b of 
the total energy is transferred by the Poynting flux according 
to Equation ( |20| |. the magnetic field in the comoving frame is 



B 



2eBL 



8.2- 10^ 



L 



1/2 
52 



fiir2 



G. 



(28) 



Because of strong self-absorption, most of photons are emit- 
ted at high cyclotron harmonics, i.e. at E ^ Eb where 



Ef 



heB 
2'Kmr.c 



0.95^Ii^keV 



(29) 



is the cyclotron energy in the comovi ng fra me. 

Making use of equations ( fT2] ) and ( IA21b in Appendix lA. 31 
one can compare the amount of the emitted photons with that 



in the Planck spectrum. 



1 



BB 



r ^" 
— tt 

1.6 rricC 



1.3 



.9/10 ^9/5 
'-B ^52 



isott e 



-2/5 



Eb 



1/10 13/5pl8/5-p 
'^RR 'in ^1 ^ max, 3 



9/5 



(30) 



-BB '10 ^1 

where in the last step we have used Equations (|2|, ([8]) and 
( |29] l to express 6, tt and in terms of outflow parameters. 
One sees that Comptonization of the thermal cyclotron radi- 
ation could provide complete thermalization if the following 
condition is satisfied 



Fr. 



13/18 
10 



< 11 



^1/4 r 1/2 

l/36p5/18 ■ 
^BB ^ max, 3 



(31) 



Compatibility with condition (|5j again restricts the flow 
Lorentz factor to mildly relativistic values in the thermaUza- 
tion region. 



19/31 9/62^1/2 
p < r 4 ^BB -^52 
— ^13/31 p5/31 

^rad max, 3 



(32) 



The constraint dsTl i on radii and Lorentz factors of the ther- 
malization zone is shown in Figure|4](solid lines), along with 
the constraint (|5]l from observed peak energies (dashed lines). 
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Figure 4. Constraints on r and F from the requirement of efficient thermal- 
ization by cv'^lotron emission (solid lines), and from observed peak energy 
(Equation (s), dashed lines). The blue, red and green lines correspond to lu- 
minosities L = 10^^, L = 10^^ and L = 10^"^ erg s~^, respectively. Other 
parameters: e^ad = £B = ebb = 0.5, Fmax = 10'^. For each luminos- 
ity the efficient thermalization region is below the corresponding solid fines. 
The crossing point of the sofid and dashed lines defines the maximal Lorentz 
factor at which the thermalization condition is compatible with condition (S)- 

If condition ( [3T1 ) is not satisfied and complete thermaliza- 
tion is not achieved, the position of the spectral peak can be 
estimated by equating the total dissipated luminosity to injec- 
tion rate of cyclotron photons times the energy they attain by 
saturated ComptonizatiorQ 



EBBi = ^r^ArSfeBTF^^cycl, 



(33) 



^ Henceforth we use ebb to denote the fraction of total energy dissipated 
in the photon-production zone regai'dless of whether a blackbody spectrum 
is attained or not, while e^.^^ still refers to the radiative efficiency as seen by 
the observer. 
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where Ar is range over which the dissipation operates and the 
Lorentz invariant A/ is given by equation (IA21l i. We find 

1/2 p3/2 -p5/18 

^ a-n I rj^ ^lad cnn "'^'^ '"lO ^ 1 ^ max,3 
Spk « 6r fceT = 590 —prrr^ ^77^ 



52 



-BB 



keV. 
(34) 



If we take Ar cx r, we obtain Epk oc r^/^, i.e. the peak en- 
ergy of the spectrum increases if the dissipation zone is moved 
beyond the thermalization radius. 
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Figure 5. The .spectral peak energy as a function of the dissipation radius 
for complete and incomplete themialization of cyclotron radiation. The dis- 
sipation operates in the range [rmax/2, rmax]. The solid lines are the nu- 
merical results, the dotted lines correspond to the blackbody result (4), with 
^rad = Ebb = 1/2. The dashed li nes c on'espond to the analytic result for 
incomplete thermalization. Equation j34t . 

The behaviour of the spectral peak as a function of the dissi- 
pation radius is shown in Figure[5](solid lines). At sufficiently 
small radii, thermalization is complete and the peak energy 
follows the blackbody law (dashed lines). It has a minimum at 
the radius where thermalization just begins to fail, after which 
the spectral peak starts increasing with radius, in agreement 
with the analytical result (|34] | (dotted line). 

Figure ^ depicts the spectra at the end of the dissipation 
episode, for different dissipation radii. For comparison, the 
corresponding thermal spectra of the same luminosities are 
also shown. As expected, at small radii the numerical result 
is practically indistinguishable from a Planck spectrum. At 
a larger radius where thermalization fails, two distinct com- 
ponents appear in the spectrum (green dashed line): the self- 
absorbed cyclotron emission (below ~ 1 keV), i.e. Rayleigh- 
Jeans with brightness temperature equal to the electron tem- 
perature, and a Comptonized spectrum at higher energies, 
with Le oc at i? <C fceT and a Wien bump near ksT 
as expected in a saturated regime. Note that as these spec- 
tra are produced deep in the flow where tt ^ 1, their shape 
has no direct bearing on the observed emission released at the 
Thomson photosphere. However, if the flow is radiatively ef- 
ficient at tt ~ 1, their peak energy is reasonably close to the 
final value. 

8. THERMALIZATION BY COMPTONIZING SYNCHROTRON 
RADIATION 

Synchrotron emissio n is a copious source of photons. 
[Thompson et al.l ( l2007h suggested that thermalization occurs 
via Comptonization of synchrotron photons when the jet 
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Figure 6. Spectra at the end of the dissipation episode, taking place in the 
range [rmax/2, rmax]. The solid (red), short-dashed (blue) and long-dashed 
(green) hnes coiTespond to different dissipation radii. The dotted hnes ai'e 
the con'esponding Planck spectra for the same luminosities and radii. The 
vertical dashed line corresponds to the average observed peak energy 210 
keV for Lrad,52 = erad-^52 = 1/2. 

breaks out of the Wolf-Rayet core at the distance of about 
10^° cm from the source. With this and some other assump- 
tions, they have managed to reproduce the observed relation 
between the peak energy and isotropic luminosity (Equation 
([1])). When analyzing the thermalization process, they as- 
sumed that all synchrotron photons emitted above the syn- 
chrotron self-absorption frequency are redistributed to i? 
SfceT; then even a small population of relativistic electrons 
in the presence of a moderate magnetic field suffices for ther- 
malization. However, most of synchrotron photons are emit- 
ted near the self-absorption frequency where the brightness 
temperature is relativistic (it is equal to the energy of emitting 
electrons) whereas thermal electrons have a non-relativistic 
temperature. Hence the energy could not be transferred from 
the "cold" electrons to the hot radiation. On the contrary, the 
high brightness temperature photons lose energy (via induced 
scattering, see AppendixlAt. Only photons emitted well above 
the self-absorption frequency, where the brightness tempera- 
ture is less than the temperature of thermal electrons, could 
be redistributed by the Comptonization to the thermal peak. 
Taking account of this, let us address the thermalization of 
synchrotron radiation at GRB conditions. 

The synchrotron emissivity of a relativistic power-law dis- 
tribution A^(7) — Nq'^^'p of electrons is approximately given 
by 



jsynch(-E) — jo 



E 



-(p+5)2 



TT 
— ( 

4 



Eb 



Eb 



-(p+5)/2 



(35) 



Assuming that electrons are injected with typical Lorentz fac- 
tor 7acc with power (3inj (per unit volume, in mcC^ units), 
we can determine the electron distribution function from the 
equation for electron cooling 



d_ 

dj 



[7A^(7)] + ^%-7acc) = 0, 



where 



4crT7^([/s + C/r 



ad; 



nicC 



(36) 



(37) 
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where w&Ub and f7rad are the magnetic and radiation energy 
densities, respectively. The solution is 



N{-f) = Nor 



3tocC( 



(38) 



Introducing einj as the fraction of total energy provided to ac- 
celerated particles over a dynamical time, we can write 



(39) 



Using this in Equation (l38l l and substituting No back to the 
emissivity expression i35[ . we obtain 



Jo 



Stt _i cF 
16 f 



(40) 



r 7acc(£B + £bb) V 

The distribution ( |38] | also yields p = 2 in Equation ( |35] |. 

The transition energy Eq ab ove w hich photons are upscat- 
tered is found from condition ( IA26b in Appendix lA.4l substi- 
tuting jsynch from Equations (l35b and ( |40| ). 



5/2 



32 TTloC Eb 

81 cttA^ (fcBr)2 



Jo- 



(41) 



With account of Equation (|A27| i the thermalization condition 
(fTTT ) is now written as 



-mj 



Tacc (e_B + £bb) 



-,4/5 / p X 6/5 



-13/5 > 

(42) 

Substituting 6*, tt and E's from Equations (|2|l, dSJ and (|29] l, 
we have finally 



0.8 



^1/8 p3/8 pl/4 
/acc,2 ' 10 1 



3/16 
52 



max, 3 ^ £ 



3/4 
B 



13/16 



-BB 



es + £bb 



(43) 



Note that the factor £i,ij /{eb +£bb ) in Equation (l43T l is always 
smaller than unity because £inj < ebb, as £bb includes the to- 
tal radiation energy while Sinj corresponds to the non-thermal 
dissipation channel only. Taking into account a very weak 
dependence on L, r, T and Fmax, this condition in fact con- 
strains the Lorentz factor of accelerated electrons from above. 

The arguments leading to the above constraint are valid 
only if Eq is below the synchrotron energy for electrons with 
7 = 7acc, i-e. if 



Eo/Eb < 0.37 



2 

acc ■ 



(44) 



This results in a lower limit for 7acc, found by combining 
Equations dTTT i and ( l44l i. along with 9 = Epi^/dTnicC^ and 
ii^pk from Equation ([T]i: 



>10 e. 



X e 



-1/6^-5/12^1/3 p5/6pl/6 



-1/12 



^10 ^1 ^ max, 3 
1/6 



£b + £bb 



(45) 



If condition i45[ is not satisfied, most of the synchrotron pho- 
tons are downscattered by the induced Compton process and 
thermal Comptonization will be depleted of seed photons, re- 
sulting in high-peaked spectra. In summary, the requirement 
that i?pk be kept below a certain value restricts 7acc to a range 
of values determined by Equations (l43T l and (l45T l. 



If complete thermalization is not achieved, the position of 
the spectral peak can be estimated by equating the total dis- 
sipated luminosity to the number of upscattered synchrotron 
photons times the energy they attain by saturated Comp- 
tonization: 



IGtT 2 2 7 
EbbL^^-CT r SfceTtdynA/synch. 



(46) 



This yields 



^1/7^5/7^1/7 4/7 

m r>-r-< T m ^I'ad r-r^n 52 1 max, 3 'acc,2 

-E'pk ~ 6r kBi = 560 

£"BB ' 
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.2/7 3/7 



-mj 



£b + Ebb 



-4/7 



keV. 



(47) 



In the domain of its applicability, the theoretical scaling 
£^pk oc r~^/^ agrees well with the numerical results. One 
can find a condition for the flow parameters by comparing 
Equation ( |47] | to the observed £'pk values from Equation 



p5/4 -pl/4 
2 7acc,2 i 1 1^ max, 3 ^ £ 



1/2 ,3/4 
BB 



7-5/8 3/4 
-'^52 'lO 



-rad 



£b + £bb 



(48) 



Note that the radius appears in the denominator on the left 
hand side of ( |48] ), meaning that large radii are favoured for 
producing low-peaked spectra. On the other hand, r is still 
constrained from above by the condition of saturated Comp- 
tonization, Equation (|9]). 

The synchrotron photons are complemented by cyclotron 
emission from thermal electrons. To determine which is the 
dominant source of photons for Comptonization, it is suffi- 
cient to compare the transi tion e nergies Eq for the two pro- 
cesses, given by equations ( |A20t and dTIT i. respectively. The 
process for which Eq is larger, dominates. We find that pho- 
tons are supplied by synchrotron if 



11/7^4/7 



2.3 



rii r 



r 11/14 4/7 
-^52 7; 



max. 3 ^ 3/14 
TT/7 - 



acc, 3 —2 



£b + £bb 



-inj 



4/7 



(49) 



The behaviour of the peak energy with increasing dissi- 
pation radius is plotted in Figure [T] The results with syn- 
chrotron emission from non-thermal electrons are shown by 
solid and long-dashed lines for two different electron accel- 
eration Lorentz factors, 7acc = lO'^ (solid) and 10^ (long- 
dashed), and bulk Lorentz factors F = 10 (red) and 20 (blue). 
High 7acc results in few synchrotron photons and the radius 
where thermalization fails is determined by cyclotron emis- 
sion (solid line). As the dissipation radius is increased, the 
peak energy initially behaves as it would without non-thermal 
injection (short-dashed line). Synchrotron emission becomes 
relevant near r = 10^^ cm , where the condition ( |49] | is ful- 
filled, and the increasing number of synchrotron photons leads 
to a decrease of iJpk with dissipation radius. 

Lowering 7acc increases the number of synchrotron pho- 
tons. In the case with 7acc = 10^ synchrotron emission 
becomes the dominant source of soft photons approximately 
at the radius where complete thermalization fails. Owing to 
Equation ( |47] i the peak energy only slightly deviates from the 
blackbody case at larger radii (for which E'pk.BB oc r~^/^). 
The decrease of E'pk continues as long as the Compton pa- 
rameter is sufficiently large to be able to upscatter most of 
the photons to ~ 3kBTe, i.e. Comptonization proceeds in a 
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Figure 7. The spectral peak energy as a function of dissipation radius for 
complete and incomplete thermalization of synchrotron radiation. The dis- 
sipation operates in the range [rmax/2, rmax]- The soHd and long-dashed 
Hnes correspond to the numerical results with non-thermal particle injec- 
tion with 7acc = 10^ and 10^, respectively. The dashed green line is the 
numerical result with thermal heating only. The dotted lines correspond 
to the blackbody result (4), with Eiad = £bb = 1/2. The horizontal 
dashed line corresponds to the average observed peak energy 210 keV for 

-f'rad,52 = erad-^52 = 1/2. 

saturated regime. The fact that the peak energy continues de- 
creasing with dissipation radius even after thermalization fails 
sets the synchrotron process apart from all other mechanisms 
considered in this work. This results in the least stringent 
constraints for the photon production site, which can extend 
to r ^ 10^^ cm. Furthermore, full thermalization is no longer 
required for obtaining low-peaked spectra. Nevertheless the 
peak energy compatible with the observations could not be 
obtained at 7acc = 1000, in agreement with the condition 
(US. 




log E [MeV] 

Figure 8. Spectra at the end of dissipation, rmax, from Comptonization of 
synchrotron radiation. The short-dashed (blue), long-dashed (green), sohd 
(red) and dot-dashed (magenta) lines con'espond to different dissipation radii. 
The dotted lines are the corresponding Planck spectra for the same luminosi- 
ties and radii. The vertical dashed line con'esponds to the average observed 
peak energy 210 keV for Lrad,52 = eradi52 = 1/2. 

Figures |8] |9] and [TO] depict the spectra for different injec- 
tion energies and bulk Lorentz factors (7acc, T) — (10'^, 10), 
(10^, 10) and (10^, 20), respectively. The different spectra in 
each figure correspond to different points on the red solid and 
dashed lines in Figure In addition to lowering the peak 
energy, increasing the dissipation radius leads to progressive 
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Figure 9. Same as FigurefS] but for electron injection Lorentz factor "Yacc — 
10^. The straight dotted line corresponds to EL^ oc E^-''^ , expected in 
the domain where the spectrum is determined by the combined effect of syn- 
chrotron emission and downscattering by the induced Compton process (see 
Equation (A25)). 

softening of the low-energy spectrum due to the decreasing 
j/-parameter A broad range of low-energy slopes can be 
produced, from thermal-like to almost flat in EFe, covering 
nearly the entire observed range. Furthermore, if y is modest, 
the low-energy slopes will survive up to the Thomson photo- 
sphere and directly correspond to the observed spectra. Note, 
however that if the dissipation radius is too large so that y < 
a few, the spectrum exhibits a broad maximum and resembles 
neither the canonical Band shape nor a cutoff power-law (e.g. 
the Tniax = 6 X 10^^ cm case in Figure [TOb. 
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Figure 10. Same as Figure H but for bulk Lorentz factor T = 20. 

The observed high-energy emission has to be produced at 
even larger radii, presumably close to tt ~' 1. The compact- 
ness parameter is very large in the photon-production region 
and any radiation above iricC^ (comoving frame) would be 
absorbed by pair production. Also, as we require y > \, the 
analogous parameter for downscattering xtt > y > 1, thus 
photons above the spectral peak are downscattered before they 
escape. 

In the regime where the seed photons are provided by syn- 
chrotron, two characteristic photon energies appear in the 
spectrum below the peak. The higher of those, Eq, is the en- 
ergy above which the produced photons are up- rather than 
downscattered. The lower one is the energy below which the 
emitted synchrotron photons are reabsorbed before they are 
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significantly downscattered by induced scattering. The spec- 
trum between these frequencies is determined by the com- 
bined effect of synchrotron emission and induced downscat- 
tering and is given by Equation ( IA25b : it is compatible with 
the observed low-energy part of GRB spectra. 

9. DISCUSSION 

In broad terms GRB emission models can be divided 
into two classes based on whether the radiation is pro- 
duced in optically thin regions or near and below the Thom- 
son photosphere. The most prominent dissipation mecha- 
nism in the f ormer category is the (collisionless) internal 
shock model jNaravan et alj |1992t iRees & Meszaroslll994 



iSari & Piranl 119971: iKobavashi et al.l1l997h . while numerous 
models of energy release in Foynting-dominated jets have 
also been considered (e.g.lUsov 1992; Lyutikov & Blandford 
2003- Granot et al. 201 1; Lvutikov 201 l : | Zhang & Yanll20li[ 
Granoti,2012; .McKinnev & Uzdenskv.,20121) . The latter cate- 



gory includes various m echanisms invoking dissipation of ei- 
ther the magnetic field (U hompsor 1 1 994t iDrenkhahn '2002'; 



Thompso 



3i 



2006t [Thomp son et al] |2007|; iGiannios 2008; 
Lvubarsk\ l2010h or the kinetic energy of relative motions 
within the jet, e.g. via inte rnal shoc ks jPe'er et alj 120061 ; 
iBromberg et al. 201 la; Levin sonl2012l) or collisional dissipa- 
tion (iBeloborodovi2010>;iVurm et al.ll20llh . 

A significant fraction of observed prompt GRB spectra ex- 
hibit a low energy turnover that could not be produced by 
radiatively efficient optically thin emission. Optically thick 
emission models avoid this problem, the local spectrum be- 
low the photosphere can be as hard as Rayleigh- Jeans or even 
WierQ. On the other hand, optically thick models face a dif- 
ferent non-trivial issue: the jet has to be able to generate a 
sufficient number of photons, otherwise the spectrum peaks 
at energies that are too high compared with observations. Our 
main finding in this work is that this requirement places severe 
constraints on the flow parameters irrespective of the details 
of the subphotospheric dissipation mechanism. 

Deep enough below the Thomson photosphere electrons 
and photons are in equilibrium and share a common tempera- 
ture (e.g. iGianniosL2012,) . the sufficient condition for which is 
that the Compton parameter y S> 1. However, even there the 
radiation does not necessarily have a Planck spectrum since 
the interaction between photons and particles is mediated by 
(single-) Compton scattering which conserves the number of 
photons. A successful thermalization relies upon the exis- 
tence of a sufficiently powerful photon source, i.e. an emis- 
sion/absorption process. 

The importance of the thermahzation issue lies in its con- 
nection to the total number of photons carried by the jet at any 
given radius and consequently the mean available energy per 
photon. Consistency of the general thermodynamic al relation 
between the radiation temperature (i.e. the spectral peak po- 
sition) and the energy density with the observed - L re- 
lation quite generally limits the thermalization region to r > 
a fewx 10^ cm if the flow is almost comp letely stopped at this 
radius, and to even larger radii otherwise ( lEichler & LevinsonI 
[2000; Thompson et al. 2007). This is much larger than the 
size of the central engine, thus relic photons from the center 
are insufficient to account for the observed spectral peak po- 
sitions, i.e. most of the observed photons have to be produced 

^ The radiation received by the observer is subject to softening due to ra- 
diative transfer effects, unless the flow remains radiation-dominated up to the 
Thomson photosphere (see e.g. ^eloborodov 2010, 201_1}. 



much further out. 

We explored systematically all radiative processes capable 
of photon production in GRB jets: bremsstrahlung, double- 
Compton scattering, cyclotron and synchrotron emission. In 
all cases efficient production of photons is confined to dense 
regions of the jet where the Thomson optical depth is still very 
large. Thermalization can be achieved either directly if the 
flow is optically thick to absorption or via saturated Comp- 
tonization in which case the emission/absorption process acts 
as a source of soft photons. The high density requirement sets 
an upper limit to the thermalization radius and, most impor- 
tantly, to the bulk Lorentz factor at the thermalization zone, 
which is constrained to be at most F ^ 10 at radii r ^ 10^"- 
10^^ cm. This rules out the possibility that sufficient amount 
of photons could be produced in a freely expanding flow, 
which would have acquired a substantially larger Lorentz fac- 
tor by those radii. 

Smooth unmagnetizedjets yield the strongest constraints on 
the photon production and/or thermalization conditions. For 
example, thermalization by bremsstrahlung requires the jet to 
be at most mildly relativistic at radii r ~ 10^" cm. Under 
typical conditions within GRB outflows, photon production 
by bremsstrahlung is g enerally superseded by double Comp- 
ton scattering (see also lBeloborodovll20 1 21) . However, strong 
constraints for successful thermalization are obtained also for 
this process: the Lorentz factor is limited to F < a few at the 
optimal radius r ^ a fewx 10^*^ cm and to even lower values 
at smaller or larger radii. 

Strong local density enhancements/clumping can make 
bremsstrahlung a plausible thermalization mechanism, ow- 
ing to the quadratic dependence of the rates on the density. 
However this requires regions of the flow to be compressed 
to densities several orders of magnitude higher than the aver- 
age. The maximal conceivable compression can be expected 
in current sheets in a Poynting dominated flow with an alter- 
nating magnetic field configuration, which can yield ~ 10^ 
for the ratio of densities in- and outside the current sheets. 
We find that in even in this case, consistency with the ob- 
served peak energies can be achieved only if F < a fewx 10 
at r ^ 10^^ cm. 

In unmagnetized flows bremsstrahlung and double Comp- 
ton are the only plausible photon production mechanisms. In 
magnetized jets cyclotron and synchrotron emission can also 
serve as photon sources. Thermalization by these processes 
relies on saturated Comptonization to bring the emitted pho- 
tons to the thermal peak. This requires that the Compton pa- 
rameter y ^ 1, which restricts the Lorentz factor in the ther- 
malization zone to at most ~ 20. The constraints for generat- 
ing sufficient amount of photons by cyclotron emission are no 
less restrictive than those of bremsstrahlung or double Comp- 
ton scattering. Additionally, in this case one also requires a 
strong magnetic field. In the optimal case of equipartition be- 
tween magnetic and thermal energy, one should have F < a 
few near r ^ 10^° cm to account for the observed spectral 
peaks. 

Finally, the observed photons may be produced by syn- 
chrotron emission provided that a fraction of electrons can 
be accelerated to relativistic energies. Although most of the 
emitted synchrotron photons are downscattered by induced 
Compton owing to their high brightness temperature and sub- 
sequently reabsorbed, the number of surviving/upscattered 
photons may be sufficient to account for observations. In fact, 
the constraints for the dissipation radius and Lorentz factor for 
synchrotron, F < 10-20 at r ^ 10^^-10^^ cm, are the least 
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restrictive among all processes considered. Most importantly, 
in this case a complete thermalization is not necessary for ob- 
taining low-peaked spectra. On the other hand, the mecha- 
nism requires efficient channeling of power to acceleration of 
electrons, within a narrow range of Lorentz factors, between 
7acc ~ 10 and a few hundred, as well as a strong magnetic 
field. 

9.1. Possible thermalization models 

We turn now to examine several astrophysical processes 
that have been suggested as sources of dissipation near the 
photosphere and examine whether the needed conditions for 
photon production can be achieved in these scenarios. 

The radius where a highly efficient kinetic to internal en- 
ergy conversion has to take place is in the vicinity of the 
progenitor surface. Indeed, hydrodynamical simulations of 
the jet-star interaction indicate that GRB outflows are far 
from simple passively expanding fireballs. Instead, the out- 
flow experiences strong heating via internal recollimation 
shocks arising from interaction with the star as well as the 
hot cocoon surroundin g the j et jLazzati & B egelmad 120051 : 
iBromberg & LevinsonI 120071: iBromberg et'ail 1201 Ibl) . The 
latter is built from shocked material as the jet burrows its way 
through the progenitor and can affect the outfl ow out to sig 



nificantly larger rad i i than the progenitor radius ( Lazzati et al 
120091: iMizuta et al.l 1201 lb . It was shown by Lazzati et al 



( 120091) that heating by recollimation shocks can maintain a 
high radiative efficiency all the way to the Thomson photo- 
sphere. In addition to heating the flow these shocks slow it 
down and the bulk Lorentz factor drops from a few hundreds 
to a few tens. However, the detailed properties of the radiation 
field depend on the heating history of the jet and subphoto- 
spheric dissipation doesn't automatically give ri se to a Planck 
spectrum. In the particular simulation shown by lLazzati et alJ 



spec 

m 



the recollimation shocks take place at r > 10 cm 
and the Lorentz factor is of the order of a few tens. Full ther- 
malization is not successful under these conditions and the 
spectrum is expected to peak at energies substantially higher 
than those typically observed. Note also that the collimation 
shock is expected to be radiation mediated (Bromberg et aP 
1201 lah and shares the problems of these shocks discussed be- 
low. 

An altern ative possibility for effi cient dissipation has been 
proposed bv lGhisellini et alj (120071) . In this model the outflow 
is shocked by encountering a small amount of slow material 
near the progenitor surface. For example, shells in an inter- 
mittent outflow can collide with the progenitor or cocoon ma- 
terial that has partially refilled the funnel excavated earlier by 
the jet. A mass fraction lO^'^-lO^'* of the displaced progen- 
itor matter is enough to have significant dynamical effects on 
the jet, resulting in substantial conversion of kinetic energy to 
heat and in slowing down the jet to F 10. Clearly, if such 
conditions can be achieved this can plausibly result in com- 
plete thermalization of the dissipated energy. However, the 
model is not sufficiently determined to check its consistency. 
An interesting point is that short GRBs show similar spec- 
tral features to long ones, albeit with somewhat harder peak 
energies. These short GRBs most likely don't involve col- 
lapsars (jNakar 2007). Therefore, both previously discussed 
dissipation scenarios, which are based on collapsars, are not 
applicable to short GRBs and one has to invoke a different 
mechanism to produce thermalization in these bursts. 

The canonical internal shocks occurring in variable out- 
flows can also give rise to subphotospheric heating. Such 



shocks are expected to be mediated by Compton scattering 
dLevinson & BromberdllOOSl: iBudnik etalllloTol: iKatz et all 
120101) . rather than plasma instabilities that operate in col- 
lisionless shocks. The downstream pressure in radiation- 
mediated shocks in GRBs is provided by photons advected 
from the upstream, the photon produc tion in the shock transi- 
tion layer is neg ligible in comparison (IBromberg et al.ll201 lat 
lLevinsonll2012l) . The photon number density in the imme- 
diate downstream is substantially smaller than in the Planck 
spectrum, which is compensated by a higher temperature than 
in the blackbody limit to support the required radiation pres- 
sure. Thus the model is also subject to the photon production 
problem, leading to difficulty in p roducing low-pea ked spec- 
tra. Very recently, it was shown by iLevinsonI ( 120 1 2|) that ther- 
malization further downstream is successful if the shocks take 
place at tt > IC*, which is in general agreement with our re- 
sults. However, since internal shocks are the result of relative 
motions within the jet, it is unlikely that they can decelerate 
the flow to r ~ a few at r > 10^*^ cm necessary for obtaining 
low-peaked spectra. 

Thus, we conclude that it is very difficult to find a subpho- 
tospheric dissipation model that satisfies the conditions for 
sufficient photon production in unmagnetized jets. This re- 
flects the difficulty of slowing the jet down to mildly relativis- 
tic Lorentz factors at large enough radii. 

The situation is even more speculative in the likely case 
where the flow is magnetized to some extent. Shocks are 
unlikely to be radia tively efficient even w ith a sub-dominant 
magnetic field (e.g. iNarayan et all 1201 1*). Furthermore, par- 
ticle acceleration required for thermalization by synchrotron 
emission is unlikely to work in shocks occurring deep in 
the flow irrespective of the flow magnetization. In the 
weak field limit these shocks are radiation mediated and 
cannot accelerate particles to high energies since the shock 
thickness is much larger than any kinetic scale involved, 
which leads to particle thermaliza t ion rather than acceler- 
ation ( IBromberg & LevinsonI 120071: ILevinsonI 120121) . It is 
also well known that particle accelera tion even in moder- 
ately m agnetized shocks is not ef ficient dLangdon et al.|[T988l : 
ISironi & Spitkovsky 200^1 l20Tlh . 

It has long been speculated that the GRB jets are launched 
as Poynting dominated outflows. Our results suggest that in 
this case, a significant part of the electromagnetic energy is re- 
leased at a scale of 10^° — 10^^ cm at relatively small Lorentz 
factors, much less than those corresponding to the total con- 
version of the burst energy into the kinetic energy of the flow. 
This could be achieved only via efficient magnetic dissipa- 
tion, therefore one has to conclude that the magnetic field is 
highly inhomogeneous at small enough scale so that multi- 
ple current sheets are formed within the flow. There are two 
possibilities: a) the fields of alternating polarity are imprinted 
already at the launch site so that a striped wind is formed; b) 
the initially r egular magnetic fie ld is destroyed either by MHD 
instabil ities (iLvutikov & Blandford.2003i) or by internal col- 
lisions (lZhang& Yanll2011l) . The last mechanism could work 
only if the plasma inertia is not small, i.e. if the plasma ki- 
netic energy is a significant fraction of the magnetic energy; 
this is excluded by our upper limit on the flow Lorentz factor 
There is also evidence that MHD instabilities could hardly do 
the job; at least the kink instability, which is the most dan- 
gerous one, is not disruptive i n Poynting domina ted flows 
dMcKinnev & Blandfordir2009l: iMizuno et al.ll2012l) . There- 
fore the striped wind seems to remain the only viable option. 

The striped wind structure arises naturally if the cen- 
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tral engine is an obliquely rotating magnetar (lUsovl 119921 : 
IZhang & Meszarosll2001l:lBucciantini et al.l2007h . If the jet is 
launched by rotating, accreting black hole, the structure of the 
magnetic field is less certain; possible scen arios for the forma - 
tion of the striped wind are discussed by 'Thompson ( 119941) . 
ISpruit et alJ (1200 Ih andlMcKinney & Uzden sky'(2012|). An 
important point is that the magnetic dissipation is accompa- 
nied by the flow acceleration so that most of the energy is re- 
leased when t he flow Lorentz factor is comparable with the 
maximal one (' DrenkhahnI 120021: iDrenkhahn & Spruit! I2002t 
Eyubarsky 2010^ 

However, a specific property of the mag- 
netic dissipation mechanism is that the acceleration is rather 
slow so that a significant fraction of the total energy could 
be re leased when the flow has not been maximally a cceler- 
ated (lGiannios& Spruit! l2007t iGianniosI lIOOSl lloil) . It is 
also possible that within the progenitor star, the acceleration 
rate is additionally reduced by interaction with the surround- 
ing matter Therefore our findings do not exclude the striped 
wind scenario but place severe limits that should be taken into 
account in future research. 

9.2. Conclusions 

In summary, we find that bulk of the observed photons have 
to be produced in the range 10^" ^ ^ S ^ few x 10^^ cm, at 
modest Lorentz factors of the order T ^ 10. Since the avail- 
able energy per baryon is typically much larger, this means 
that the flow kinetic energy is small in that region, and most of 
the available energy is in internal form (magnetic or thermal). 
Photon production by processes involving only thermal elec- 
trons/pairs yield the strongest constraints, limiting the flow 
to be at most mildly relativistic in the thermalization zone. 
The only exception is bremsstrahlung emission from strongly 
clumpy flows, in which case the bulk Lorentz factor must still 
satisfy F < 30. Synchrotron emission from relativistic elec- 
trons can be a copious source of photons and relax the limit 
on the Lorentz factor to F < 20. However, it requires very 
specific conditions: a highly efficient non-thermal accelera- 
tion mechanism transferring a significant fraction of the total 
energy to electrons with Lorentz factors between 10 and a few 
hundreds as well as a magnetic field near equip artition. 

The constraints are somewhat less restrictive if the radi- 
ation experiences significant adiabatic losses on the way to 
the Thomson photosphere, which might be the case if the 
flow reaches the coasting stage before transparency. Adia- 
batic cooling lowers the mean photon energy, however this 
happens at the expense of the radiative efficiency. Observa- 
tionally, comparison of the energy emitted in the prompt and 
afterglow phases indicates that a large fraction of the avail- 
able energy is typically released as prompt emissio n (e.g. 
Freedman & Wa xman 2001; Granot et al. 2006; Fan & Piran! 
20061: IZhang et alj 120071) . implying high efficiency. There 
is also evidence of a positive correlation between the radia- 
tive efficiency of the prompt phase and the burst luminos- 
m; (iLlovd-Ronning & Zhan j !2004l; !Eichler & Jontof-Huttei! 
|2005; Margutti et al.l !2012!) . The lower efficiency of less 
luminous lower-peaked bursts may alleviate the photon- 
production problem which is otherwise more serious for low- 
luminosity bursts. 

Substantial decrease of the jet opening angle between the 



thermalization region and the Thomson photosphere can also 
relax the co nstraints on the jet p arameters at the photon pro- 
duction site ( lBeloborodov!l20 1 2|) . Most of this range lies be- 
yond the progenitor surface where the opening angle is un- 
likely to change by a large factor, thus the conclusions reached 
in this work should not be significantly affected. 

Although in this work we addressed only the problem of the 
spectral peaks, it is clear that the entire observed spectrum and 
in particular the high-energy part cannot originate from the 
thermalization region. At these radii the radiation compact- 
ness is huge and any high-energy photons would be imme- 
diately absorbed by photon-photon pair production. Further- 
more, the high Compton parameter ensures that any photons 
above the thermal peak will be downscattered even if they are 
below the pair-production threshold and can't be absorbed. 
Thus the production of photons and their redistribution into 
non-thermal spectra must be the result of either two different 
dissipation episodes or a single continuous one operating over 
a very broad range of radii. Note in particular that (thermal) 
Comptonization that should take place beyond the thermal- 
ization zone and produce the high energy tail cannot produce 
photons at energies higher than ~' TmcC^. Thus, the highest 
energy photons must be produced by yet another, non-thermal 
mechanism giving rise to a relativistic electron population, 
operating either in the prompt or the afterglow phase. 

In baryonic jets, dissipation by neutron-proton collisions 
is a promising mechanism for producing the high-energy tail 
([geloborodov 2010): the requirement that bulk of the energy 
be in internal form at r ~ 10^*^ along with the low F facil- 
itates the formation of a compound flow (counter-streaming 
proton and neutron components). Compared to a passive fire- 
ball, the saturation radius is now pushed further out relative to 
the neutron-proton decoupling radius, i.e. the flow may still 
be accelerating when neutrons and protons start migrating rel- 
ative to each other, resulting in strong collisional heating. In- 
ternal shocks can also provide additional energy release near 
the Thomson photosphere. 

In Poynting-dominated jets, gradual dissipation of the 
"free" magnetic energy associated with small-scale field re- 
versals can be responsible both for photon production as 
well as distorting the spectrum into a non-thermal shape. 
Heating in current sheets maintains a thermal spectrum un- 
til the bremsstrahlung photosphere, while continuing dissi- 
pation at larger radii can give rise to the observed extended 
high- energy component by Co i nptonizing the b lackbody pho- 
tons (!Giannios & Spruit!!2007t !Gianniosl!2008l) . Whether the 
mechanism can reproduce the observed spectral peaks de- 
pends on the ability of the jet to keep its Lorentz factor low 
enough to sufficiently large radii. 
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APPENDIX 

A. COMPTONIZATION OF SOFT PHOTONS 

In this Appendix, we address the evolution of the radiation spectrum in the comoving frame. Well inside the scattering pho- 
tosphere, the radiation is isotropic in the comoving frame therefore the evolution of the photon spectrum is governed by the 
Kompaneets equation 

where n is the photons occupation number, j{E) and a{E) are the photon emissivity and opacity, respectively. The last term 
accounts for adiabatic coohng; in a spherical outflow with the Lorentz factor F, 

(V-v)-^. (A2) 

In the Comptonization operator, the first term in the brackets describes redistribution of soft, E < k^T, photons towards E ^ 
SfceT, the second describes the recoil effect negligible for soft photons and the third term is responsible for the induced scattering. 

At the condition (|6]l, the Bose-Einstein spectrum with the peak at E ^ S/sbT is formed. At the sam e co ndition, one can neglect 
the adiabatic losses (the last term in Equation (lAlb ). Moreover, instead of solving the full Equation (lAlb ) one can now consider 
only the total photon balance, 

dM 

where 



A/", (A3) 



= TvW / E'^ndE (A4) 
[chy J 

is the total number of photons in the Bose-Einstein spectrum, M is the rate of injection of photons into the thermal peak. The 
latter is determined by the processes at small photon energies, E ^ k-QT; by virtue of the condition (|6]l the Kompaneets equation 
is reduced in this range to 



which means that the injected photons are steadily redistributed in energies by the Comptonization process. The photon injection 
rate could be found from the solution to this equation as the photon flux in the energy space towards higher energies but one can 
find simple estimates just by comparing different terms in this equation. 

Note that the induced scattering (the term) dominates only if the radiation brightness temperature, Tb ~ En/k-Q, exceeds 
the electron temperature, T; in this case the photons are redistributed towards smaller energies until they are eventually absorbed. 
Such a situation arises if the synchrotron emission serves as a photon source. If the emission/absorption processes are only due to 
the thermal electrons, the thermodynamical condition, j = {k-oT / E)a, implies that the radiation brightness temperature does not 
exceed that of electrons. In this case the induced scattering becomes only marginally important when the spectrum approaches 
that of Rayleigh-Jeans so that for the estimates, one can neglect the induced scattering at all and write Equation (IA5I) as 



TOcC E^ oE oE \ E ) 

Since the opacity increases with decreasing photon energy, the right-hand side of Equation ( IA6b dominates at small enough E; 
the Rayleigh-Jeans spectrum is established in this band. At larger photon energies, the Compton redistribution rate exceeds the 
absorption rate. The boundary energy, Eq, is found by comparing the left-hand side and the right-hand side of equation (IA6b as 

kvT 

^-^otN = a{EQ). (Al) 

This relation could be also found from the following considerations. In order to avoid reabsorption, an emitted photon has to 
experience significant upscattering before it is absorbed. Thus if the Compton parameter y (see equation (|6]l) associated with the 
absorption time-scale, tabs = 1 /ct{E), is greater than unity, the photon has enough time to gain a significant energy and, since a 
decreases with E, such a photon avoids reabso rptio n altogether and eventually reaches the thermal peak. In the opposite case, the 
emitted photons are absorbed. The condition ( IA7I ) corresponds to the boundary, y = \. 
Now let us consider specific mechanisms of the photon production. 

A.l. Bremsstrahlung 

At GRB conditions, bremsstrahlung can serve as a photon source only if plasma is heavily loaded by elec tron positron pairs 
(see SectionO therefore we use the bremsstrahlung opacity for the pair plasma, which can be written as (e.g.. lPozdnyakov et alJ 
[l983;.Haug„1985l) 
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or in the Rayleigh- Jeans regime 



%(^)= 31/2^3/2 [—] ^c-TcN^N., (A9) 



where Ac is the Compton wavelength and N± are the electron/positron number densities. For the electron-proton plasma, the 
opacity is just 

where A^p is the proton density. The Gaunt factor is approximated by (e.g. lRvbicki & Lightman|[l979l) 

gff{E,T)^ — In . (All) 

From Equations iAli and iA9i we find the equation for the boundary energy Eq assuming = iV+ 

UbtJ 3i/2 8 7r3/2 [m^c^J ^--^c- (A12) 

All the pho tons emitted at E > Ep are redistributed to the thermal peak so the photon injection rate is obtained as 
dlUarionov & Siuniaev.,1975;,Pozdnvakov et a l. 1983) 

M^=-% r^\^jiE)dE = 4n(J^Ya^{E = k^T)\n'^^^^. (A13) 

[chy Jeo \meC^XcJ " Eq 

A.2. Double Compton scattering 

If the radiation energy density is large, double Compton scattering becomes an important source of soft photons for Comp- 
tonization. The effective opacity for a Maxwellian electrons interacting with a Wien distribution of photons can be written at 
E/ksT < 1 as (,Thorne,. 1981: .Lightmaa 198k .Pozdny akov et al. 1983; Svensson 1984) 

aME) = ^(^) e-^g^c{0)\laTcNN, (A14) 



TT^ \ E ^ 

where — k-QT/rricC^ <^1, N and J\f are the electron and photon number densities, respectively, and 

.900(6*) = (1 + 13.9161 + 11.056(2 + 19.926'3)-^ (A15) 
is a fitting formula to the exact numerical result. For the Planck spectrum, one gets 

38.4af AbT" ^ 



"DC(S) = ^ 9BciO)aTcN. (A16) 

The calculation of the photon injection rate is similar to that for bremsstrahlung. The produced photons are upscattered by 
single Compton scattering before reabsorption above the energy Eq, found from Equation jAlj as 

''°^'=^%Dc(e). (A17) 



knT 



Now the photon injection rate is found as 



A/dc = 8^ " , aT,c{E = ksT) In (A18) 

A3. The nnal cyclotron emission 

The cyclotron radiation spectrum at high harmonics, E ^ Eb, could be considered as continuous because the cyclotron line 
width becomes comparable with the line separation. Since the electron temperatures relevant to our problem are of the order 
of 10 keV, we can use the Trubnikov approximation obtained for parameters appropriate for thermonuclear fusion dTrubnikovl 
fT958h . 

^aai'a^Nc'^e'^ (^^J . (A19) 

In lTrubnikovl d 1 95 8h the constants take the values a = 3 (1207r)2/2, q ~ 3 and s = 6, which provides a good fit at low harmonics. 
At the energies relevant to the present problem, E/Eb ~ 10, the exact opacity already exhibits a significant downward curvature 
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and the Trubnikov approximation gives a poor fit. We find that in the relevant parameter region a = 10^", 9 = 4 and s = 10 
provide a sufficiently accurate fit for our analytical estimates. 

The cyclotron emissivity and absorption coe ffici ent decrease very rapidly with frequency, thus almost all the radiation is 
emitted near the energy Eq defined by equation ( IA7b . where the upscattering and reabsoprtion rates are equal. Making use of the 
approximation ( IA19I ), one finds 

At E < Eq, the radiation spectrum is Rayleigh- Jeans. At i? > Eq, the cyclotron emission and absorption could be neglected 
and the spectrum is formed by Comptonization. The injection rate of photons into Comptonization can be found from the 
Kompaneets equation (lAll i by considering the photon flux in energy space at Eq, 



Here we substituted the photon occupation number as 



n-'S^, <A22, 

because n cx E^^ is a solution to Equation (IA6I 1 when the right-hand side is negligibly small; the normalization is chosen such 
that the spectrum goes to the Rayleigh- Jeans si E = Eq. Note that this solution describes a constant photon flux in the energy 
space from E = EqXo higher energies. 

A.4. Synchrotron emission 
The synchrotron emissivity in GRB jets could be presented as (see discussion in Section[8]l 

j(i^)=Jo(^)'''. (A23) 

The synchrotron absorption becomes important only at a very low photon energies, where the radiation brightness temperature is 
comparable with the energy of emitting relativistic electrons. At the energies larger than the self-absorption energy, the brightness 
temperature decreases but while it exceeds the temperature of thermal electrons, the Comptonization is determined by the induced 
scattering therefore the Kompaneets equation ( IA5l l is reduced to 

o^N 1 d 



E^n^^j{E). (A24) 



The solution has a form dSvunvaevI 19711) 



Substituting this expression back into the Comptonization operator, one ensures that the photon flux in the energy space is directed 
towards lower energies where these photons are eventually absorbed. 

The Compton upscattering begins when the brightness temperature decreases down to the electron temperature. The boundary 
energy, Eq, may be found from the condition that the two terms in the Comptonization operator (the left-hand side of Equation 
( IA5I )) becomes equal. For the spectrum ( IA25I ), this condition can be written as 

|fl^<^^.„^.), (A26, 

61 nicC Eq 

Note that if the emissivity were provided by the thermal electron population, j = {U-qT / E)a, the Equation ( IA26b would be 
reduced, to within a factor of about unity, to Equation ( IA7b used when only thermal electrons were considered. All photons 
emitted si E > Eq are redistributed, by virtue of the condition towards the thermal peak therefore the photon injection rate 
is found as 

,A27) 

{chf Je„ {chf 
B. NUMERICAL METHOD 

The Kompaneets equation (lAlb is coupled to the corresponding kinetic equation governing the evolution of the electron distri- 
bution, which have to be solved simultaneously for a self-consistent solution. However, in the present case the problem simplifies 
owing to the fact that thermalization can only take place well below the Thomson photosphere where bulk of the electrons are 
MaxwelUan. Furthermore, the (thermal) lepton population has negligible heat capacity and is therefore kept in a quasi-equiUbrium 



16 



VURM, LYUBARSKY, & PIRAN 



temperature determined by the combined effect of dissipation heating and radiative cooling. Thus solving a separate kinetic equa- 
tion for thermal electrons is not necessary. If the dominant cooling mechanism is Compton scattering, the electron temperature is 
found by equating the volume heating rate of electrons, Ph, with the energy transfer rate determined by the Kompaneets kernel. 
This yields 

^ ^ A3 Ph/(87rmec3gTiV) + Jx^n{l + n)dx ^^^^ 
A J ndx ' 

where x = E/mc<? is the dimensionless photon energy. The volume heating rate Ph is determined by assuming that a fraction 
Ebb of the total available energy is dissipated between rmin and r^^x = '^min + Ar from the central source. For constant power 
per logarithmic radius interval, this gives 



At each timestep. Equations dBlb and dTji together with the radiation field given by the Kompaneets equation dAll i determine the 
electron temperature and density and thus the entire thermal distribution. 

If we allow a fraction of leptons to be accelerated to relativistic energies (see Section [8] on thermalization by synchrotron 
emission), we need a kinetic equation describing their cooling. 



dt \ 2 (?7 

where N{j) is the electron distribution function, 7 and describe the cooling/heating and diffusion in energy space due to 
synchrotron emission/absorption, as well as Compton scattering in the Thomson regime. The source and sink terms jc and 
account for Compton scatterings in the Klein-Nishina regime, the former also includes the injection term accounting for 

accelerated electrons. 

To s olve the kinetic equations (lAlb and ( IB3b we use a numerical code developed bv lVurm & PoutanenI (120091) and lVurm et al.l 
(I2OI lb for simulating radiative transfer in relativistic flows. For the present problem the code has been modified in three main 
respects: first, we have included induced Compton scattering which was not present in the original version. This is crucial for 
obtaining true blackbody spectra (instead of Wien). Furthermore, induced scattering dominates ordinary scattering in cases when 
the radiation brightness temperature exceeds 6, which can be relevant if a suprathermal population of emitting leptons exists. 
Secondly, we solve the equation for the angle-averaged photon distribution instead of the full radiative transfer equation. This 
is justified by the huge Thomson optical depth in the thermalization zone that keeps the radiation field very nearly isotropic in 
the comoving frame. Finally, the electron equation is solved only for the non-thermal population (when applicable). This avoids 
numerical difficulties at low energies that are associated with keeping track of energy conservation in near-equilibrium situations, 
where large cancelling terms can lead to errors as well as numerical instabilities. These last two simplifications allow us to use 
the familiar Kompaneets kernel in the photon equation, as well as to shorten significantly the computation time. 

C. CONSTRAINT ON THERMALIZATION LOCATION IN A COLLIMATING FLOW 

Consider a jet with a total (kineticH-internal, beamed) luminosity Lt and beaming solid angle Ail(r) that may vary with radius. 
The isotropic equivalent luminosity is now also a function of radius and is given by L{r) — Lt/Af2(r). The peak energy of a 
blackbody spectrum is given by 

E,^ = 2m. f^'-^ ( keV. (CI) 
V '"12 EBB V An(r) / 

The Yonetoku relation now becomes 

T \ 1/2 

£rad -'^1,52 



where r,t is the radius where the radiation decouples from the flow and gives rise to the observed emission. Equating the two 
expressions for Ep\^, we find the constraint for the thermalization location 



r p3/2 t,52 [^^(^)]l/2 ^3/2 



cm. (C3) 



W hether the photons from the central source are sufficient to account for observed Spk-s is determined by writing Equation 
iC3\ for the base of the flow, i.e. setting r = tq and F = Fq ~ 1. Assuming a quasi-isotropic flow at the base, AQ{ro) = 2tt, we 
find the required opening angle of the jet at the photosphere 

0*~3x 10-3 fek ^0,7 ^52'^), (C4) 
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where we have used Aft — n6^. Note that thi s als o requires F > l/9i, since radiation cannot be beamed into a cone narrower 
than 1 /F. The Oi, values imp osed bv Equation (IC4[l are at least an order of magnitude below the typical opening angles inferred 
from afterglow observations (iLiang et al.ll2008t ICenko et al.ll2010l) strongly suggesting that the photons from the central engine 
are insufficient to account for the observed spectral peaks and that most of the observed photons must be produced in the jet at 
larger radii. Conversely, for a given opening angle 6^, the radius at the base of the flow has to satisfy 

ro = 3 X 10^ ^ L^^'/^ir^) cm, (C5) 

^BB 

which is significantly larger than the size of the central engine. Furthermore, the minimal value given by Equation ( IC5I ) assumes 
that the flow is quasi-spherical at tq. This is implausible at radii indicated by Equation ( IC5I ). thus even if the flow is almost 
completely stopped at vq, a more realistic minimal radius for this to happen is tq > 10^ cm. Even more extreme constraints are 
obtained for magnetized jets, in which case the factor ebb relating the radiation luminosity to the total luminosity at the base is 
much smaller than unity. 

It is reasonable to expect that most of the flow collimation takes place at relatively small radii in a dense environment providing 
strong confinement for lateral expansion. At radii larger than the progenitor size, r > 10^" cm, the angular factor in Equation 
iC3\ should not affect the constraint on r/F by more than a factor of a few. 



Abramowicz. M. A., Novikov, I. D., & Paczynski, B. 1991, ApJ, 369, 175 
Beloborodov, A. M. 2010, MNRAS, 407, 1033 
— . 2011. ApJ, 737, 68 
— . 2012, arXiv: 1207.2707 

Bromberg. O., & Levinson. A. 2007. ApJ, 671, 678 
Bromberg, O., Mikolitzky, Z., & Levinson, A. 2011a, ApJ. 733. 85 
Bromberg. O., Nakar, E., Piran, T.. & Sari, R. 2011b. ApJ, 740, 100 
Bucciantini, N., Quataert, E.. Arons, J., Metzger, B. D., & Thompson, T. A. 

2007, MNRAS, 380. 1541 
Budnik, R., Katz, B., Sagiv, A., & Waxman, E. 2010, ApJ, 725, 63 
Cenko. S. B., Frail, D. A., Harrison, R A., et al. 2010, ApJ, 711, 641 
Cohen, E., Katz, J. 1., Piran, T., et al. 1997, ApJ, 488, 330 
Crider, A., Liang, E. P., Smith, I. A., et al. 1997, ApJ, 479, L39 
Daigne, R, Bosnjak, Z., & Dubus, G. 201 1, A&A, 526, Al 10 
Drenkhahn, G. 2002. A&A, 387, 714 
Drenkhahn, G., & Spruit, H. C. 2002, A&A, 391, 1141 
Eichler, D. 1994, ApJS, 90, 877 
Eichler, D., & Jontof-Hutter, D. 2005. ApJ, 635, 1182 
Eichler, D., & Levinson, A. 2000, ApJ, 529, 146 
Pan, Y., & Piran, T. 2006, MNRAS. 369, 197 
Preedman, D. L., & Waxman, E. 2001, ApJ. 547. 922 
Ghirlanda, G., Celotti, A., & Ghisellini, G. 2003. A&A, 406, 879 
Ghirlanda. G., Nava, L., & Ghisellini, G. 2010, A&A, 511, A43 
Ghirlanda. G., Nava, L., Ghisellini, G., Celotti, A., & Firmani. C. 2009. 

A&A, 496, 585 
Ghisellini. G., & Celotti, A. 1999, ApJ, 511, L93 
Ghisellini, G., Celotti, A., Ghirlanda, G., Firmani, C, & Nava, L. 2007, 

MNRAS, 382, L72 
Giannios, D. 2006, A&A, 457, 763 
— . 2008. A&A, 480, 305 
— . 2012, MNRAS, 422, 3092 
Giannios, D., & Spruit, H. C. 2007, A&A. 469, 1 
Goldstein, A., Burgess, J. M., Preece, R. D., et al. 2012, ApJS, 199, 19 
Goodman. J. 1986, ApJ, 308, L47 
Granot, J. 2012. MNRAS, 421, 2467 

Granot, J., Komissarov, S. S., & Spitkovsky, A. 2011, MNRAS, 411, 1323 
Granot, J., Konigl, A., & Piran, T. 2006, MNRAS, 370, 1946 
Grimsrud, O. M., & Wasserman, I. 1998, MNRAS, 300, 1158 
Gruber. D., Greiner, J., von Kienlin. A., et al. 201 1, A&A, 531, A20 
Haug, E. 1985, A&A. 148. 386 

Illarionov, A. R, & Siuniaev. R. A. 1975, Soviet Ast., 18, 413 
Kaneko, Y, Preece, R. D., Briggs, M. S., et al. 2006. ApJS, 166, 298 
Katz, B., Budnik, R., & Waxman, E. 2010, ApJ, 716, 781 
Katz, J. 1. 1994, ApJ, 432, L107 
Kobayashi, S., Piran, T., & Sari. R. 1997, ApJ. 490, 92 
Langdon, A. B., Arons, J., & Max. C. E. 1988, Physical Review Letters, 61, 
779 

Lazzati. D., & Begelman. M. C. 2005, ApJ, 629, 903 

Lazzati, D., Morsony, B. J.. & Begelman, M. C. 2009. ApJ, 700, L47 

Levinson, A. 2012, arXiv: 1205.3227 



Levinson, A., & Bromberg, O. 2008, Physical Review Letters, 100. 131101 
Li, C. & Saii, R. 2008, ApJ, 677, 425 

Liang, E., Kusunose, M., Smith, I. A.. & Crider, A. 1997, ApJ, 479, L35 
Liang, E.-W., Racusin, J. L., Zhang, B., Zhang, B.-B., & Buitows, D. N. 

2008, ApJ, 675, 528 
Lightman, A. R 1981, ApJ, 244, 392 
Lloyd-Ronning, N. M., & Zhang, B. 2004, ApJ, 613, 477 
Lyubarsky, Y. 2010, ApJ, 725, L234 
Lyutikov, M. 2011, MNRAS, 411, 422 

Lyutikov, M., & Blandford, R. 2003, arXiv:astro-ph/03123?r 

Maigutti, R., Zaninoni, E., Bernardini, M. G., et al. 2012, arXi v: 1203. 10591 
McKinney, J. C, & Blandford, R. D. 2009, MNRAS, 394, L126 
McKinney, J. C, & Uzdensky, D. A. 2012, MNRAS, 419, 573 
Meszaros, P, & Rees, M. J. 2000, ApJ, 530, 292 
Mizuno, Y, Lyubarsky, Y, Nishikawa, K.-L, & Hardee, R E. 2012, 

arXiv: 1207.4949 
Mizuta, A., Nagataki, S., & Aoi, J. 201 1, ApJ, 732, 26 
Nakar, E. 2007, Phys. Rep., 442, 166 

Narayan, R., Kumar, P, & Tchekhovskoy, A. 201 1, MNRAS, 416, 2193 
Narayan, R., Paczynski, B., & Piran, T. 1992, ApJ, 395, L83 
Paczynski, B. 1986, ApJ, 308, L43 

Pe'er, A., Meszaros, R, & Rees, M. J. 2006, ApJ, 642, 995 

Pozdnyakov, L. A., Sobol, I. M., & Syunyaev, R. A. 1983, Astrophysics and 

Space Physics Reviews, 2, 189 
Preece, R. D., Briggs, M. S., Mallozzi, R. S., et al. 1998, ApJ, 506, L23 
Rees, M. J., & Meszaros. R 1994, ApJ, 430, L93 
Rees, M. J., & Meszaros. P 2005, ApJ, 628, 847 

Rybicki, G. B., & Lightman, A. P. 1979, Radiative processes in astrophysics 

(New York: Wiley-Interscience) 
Sari, R., & Piran, T. 1997, ApJ, 485, 270 
Shemi, A., & Piran, T. 1990, ApJ, 365, L55 
Sironi, L., & Spitkovsky, A. 2009, ApJ, 698, 1523 
— . 2011, ApJ, 726, 75 

Spruit, H. C, Daigne, R, & Drenkhahn, G. 2001, A&A, 369, 694 
Svensson, R. 1984, MNRAS, 209, 175 
Syunyaev, R. A. 1971, Soviet Ast., 15, 190 
Thompson, C. 1994, MNRAS, 270, 480 
— . 2006, ApJ, 651,333 

Thompson, C, Meszaros, P, & Rees, M. J. 2007, ApJ, 666, 1012 
Thome, K. S. 1981, MNRAS, 194, 439 
Trubnikov, B. A. 1958, Soviet Physics Doklady, 3, 136 
Usov, V. V. 1992, Nature, 357, 472 

Vurm, 1., Beloborodov, A. M., & Poutanen, J. 201 1, ApJ, 738, 77 
Vurm, 1., & Poutanen, J. 2009, ApJ, 698, 293 

Yonetoku, D., Murakami, T., Nakamura, T, et al. 2004, ApJ, 609, 935 
Zhang, B., Liang, E., Page, K. L., et al. 2007, ApJ, 655, 989 
Zhang, B., & Meszaros, P 2001, ApJ, 552, L35 
Zhang, B., & Yan, H. 201 1, ApJ, 726, 90 



